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A review of the status of the Monte Carlo event generator PHOKHARA, developed for experiments using the 
radiative return method. The four-pion production in electron-positron annihilation and in r-lepton decays and 
the narrow resonances studies are described. 



1. Introduction 

The basis of the radiative return method, pro- 
posed in pQ, is an observation that one can 
extract the hadronic cross section (e + e~ — > 
hadrons) from the measurement of the cross 
section of the reaction e + e~ — > hadrons + 
photons, where the photons are emitted from 
the initial leptons. This is possible due to the 
factorization da(e + e~ — > hadrons + Jisr) = 
H(Q 2 ,8- l )da(e + e~ — ► hadrons)(s — Q 2 ), where 
the function H is fully calculable within QED and 
Q 2 is the invariant mass of the hadronic system. 

The radiative return method has been used by 
meson factories DAPHNE, BaBar and BELLE 
and allows for the measurement of the hadronic 
cross section from the nominal energy of these 
experiments down to the production threshold. 

The traditional way of measuring of the 
hadronic cross section via the energy scan needs 
dedicated experiments, while using the radiative 
return one can profit from the existing meson fac- 
tories. The smaller cross section of the radiative 
process (by a factor of a/ir as compared to the 
process without photons emission) has to be com- 
pensated by higher luminosities of factories. 

To obtain the hadronic cross section using the 
radiative return method in a realistic experimen- 
tal situation, one needs a Monte Carlo event gen- 
erator of the measured process. To provide such 
a tool to the experimental groups PHOKHARA 
event generator was constructed. The construc- 
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tion of the PHOKHARA event generator started 
from the EVA [213] generator, where structure 
function method was used to model multi-photon 
emission. The physical accuracy of the program 
was however far from the demanding experimen- 
tal accuracy. To meet this ever growing demands, 
the event generator PHOKHARA [3] was con- 
structed. It is based on the complete calcula- 
tion of radiative corrections to the ncxt-to-leading 
order for the ISR emission and relevant next-to- 
leading order corrections to the final state emis- 
sion. In this paper the latest updates of the four- 
pion channels [5] are briefly outlined (Section 2) 
and preliminary results for the narrow resonances 
implementation are presented (Section 3). 

2. The four-pion production in r decays 
and e + e annihilation 

The four-pion production in e + e~ annihilation 
was implemented in the generator EVA [5] a long 
time ago and recently it was reanalyzed in [5] . 

There are altogether four different channels ac- 
cessible for 47T production: 

e + e~ — * 27T + 27r~, 

e + e~ — * 27r°7T + 7r~ (a), 

r" i/ r 27r"7r+7r (b), 

t~ — > v T 3ir ir~ . 

Assuming isospin symmetry, the amplitudes of 
either (a) [6] or (b) [7] are sufficient to determine 
all four amplitudes. 

In EVA and PHOKHARA a choice of [g was 
adopted and a function 3 ^ (symmetric under the 
interchange of p± and pi and antisymmetric under 
the interchange of p + and p~) is used to model 
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the two e + e four pion channels: 

(tT+TT-T^VmI ) - J»(P1,P2,P + , P -), (1) 

The other matrix elements: (7r^ 7r^ 7r^ | 1 0) , 
(7r-7r?7r^7r§|J-|0) and (71-f tt^tt+tt ] J~ |0) can be 
expressed as sum of J M functions with permuted 
arguments [BJ. 

For e + e~ annihilation the current J M contains 
the complete information about the hadronic 
cross section through 

J J e ; n ( J v m y d<j> n (Q;qi,...,q n ) 

= ^-{QnQv-g^vQ 2 ) R(Q 2 ) , (2) 

where R(Q 2 ) = cr(e + e~ — ► hadrons){Q 2 ) /a po i n t- 
Similarly for r decay we have 

J J^ J y* d$ n (Q; (71, ... , g n ) 

= ^-{QvQ»~9^Q 2 ) R T (Q 2 ) (3) 

where R T can obtained from the differential r de- 
cay rates. 

From the isospin relations between the matrix 
elements of the four pion hadronic current one 
obtains the relations between r decay rates and 
e + e~ annihilation cross sections 

R T {- 0) = !#(+ + --) 

JT(- - +0) = i + - -) 

+ JZ(+ - 0) . (4) 

They allow for direct tests of the isospin sym- 
metry provided all R functions were measured. 

The latest (still preliminary) accurate measure- 
ment (by means of the radiative return method) 
by BaBar [8] of the 27r°7r + 7r~ mode comple- 
mented the precise e + e~ data sets by CMD2 [9], 
SND [10] and BaBar p]. Together with the r 
data from ALEPH [H] and CI. IX) !.'. it allows 
for model independent isospin symmetry tests. 

Combining the results from BaBar [1X18] and 
using relations between r and e + e~ Eq.Q, we 
obtain predictions for the r spectral functions v 
(related directly to R T [S]). We use normalization 
of the spectral functions chosen by ALEPH [T3] . 
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Figure 1. The spectral function of the t _ — * 
3ir°Tr~v T decay mode. ALEPH [12] data versus 
predictions from BaBar [ll|8j and the model [5]. 
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Figure 2. The spectral function of the t~ — > 
27T-7r + 7r°i/ T decay mode. ALEPH [12] and CLEO 
[13] data versus predictions from BaBar [ll|8j and 
the model [5]. 
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Figure 3. Diagrams contributing to the hadronic 
current in [3] (upper) and the new contributions 
from p mesons (middle) and the omega (lower). 



As shown in Fig[T] and Fig[5] there is a good 
agreement between the spectral functions pre- 
dicted from the BaBar data and the isospin sym- 
metry and the ones measured by ALEPH and 
CLEO. One observes systematic shifts which are 
however well contained within current error bars 
and it is not possible to claim an observation of 
the isospin symmetry violation. 

Old 47r model adopted from [T3] and used in [3] 
cannot reproduce new and more accurate data. 
The update to the model has been constructed 
in [5]. Amplitudes used in the new model are 
schematically depicted in Fig[3J The upper dia- 
grams show old contributions from [3J, the mid- 
dle diagrams represent newly added contribu- 





Br("r _ -> i/ T 27r-7r+7r u ) 


pdg m 


(4.46± 0.06)% 


model 


(4.12 ± 0.21)% 


BaBar (CVC) 


(3.98 ± 0.30)% 



Br(r — > v t tt u)(n tt + tt {> )) Br(r > v T ir 37r u ) 

(1.77± 0.1)% (1.04± 0.08)% 

(1.60 ± 0.13)% (1.06 ± 0.09)% 

(1.57 ± 0.31)% | (1.02± 0.05)% 

Table 1 

Branching ratios of r decay modes. Comparison 
between model, experimental data [15] and pre- 
dictions based on BaBar data [ll|8j and isospin 
symmetry. 



tions where the p particles are treated as SU(2) 
gauge bosons and the lower diagram represent 
new omega contributions, which substitute the 
last diagram of the first line. 

There are following parameters in the model: 
external masses and widths m p >, r p i,m p n, T p », 
rripiii , Y pin , four couplings in each of the a\ , /o 
and lo parts and one coupling in p-part. The pa- 
rameters were fitted to the existing data. The fit 
is quite good, with x 2 / n d.o.f — 275/287. 

The comparison between new model and the 
data from ALEPH, CLOE and BaBar are shown 
in FigQ] and Fig[2l The upper and lower curves 
represent error bars. 

It is interesting to compare also the r branch- 
ing ratios from the PDG [15], the new model pre- 
dictions and the direct predictions from BaBar 
data using isospin symmetry. The results are 
collected in Table [1] They are in agreement 
within current error bars. One can observe about 
two sigma difference between PDG and BaBar 
data for Br(r~ — > ^ r 27r~7r + 7r°). Better preci- 
sion BaBar data, expected after the preliminary 
results [8] are published, might shed light on the 
isospin violation in the four pion e + e~ production 
and r decays. 
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3. The narrow resonances J/ip and i/j(2S) 
in PHOKHARA 

The implementation the narrow resonances 
3/tp and ^(2S) to the PHOKHARA event gen- 
erator is the latest update of this program. The 
final results will be presented soon and in this 
paper we present the preliminary results. 

Contributions from two narrow resonances 



J/ip m = 3096.916 MeV 
ip(2S) -► m = 3686.093 MeV, 

to the final states: 



T = 93.4 keV 
T = 337 keV, 



7T+7T-, fi+fjT, K+K- and K°K° 
were implemented in PHOKHARA. 

In general one has three types of contribu- 
tions to the production amplitude at the narrow 
resonance (Figf4|): one-photon continuum, one- 
photon decays and three-gluon decays. The last 
amplitude contributes only to the K + K~ and 
K°K° production. 




ggg 

Figure 4. The Feynman diagrams of e + e~ — > 
/x + /Lt _ , 7r + , 7r~, KK at charmonium resonance, 
the one-photon continuum process, the one- 
photon decays and the three-gluon decays (only 
for kaons). 



For reliable predictions of the tt + -k and KK 



production good models of electromagnetic pion 
and kaon form factors are needed. The form fac- 
tors used in the public version of PHOKHARA 
6.0 were taken from [IB]. The CLEO-c measure- 
ment |17) and results of the investigations in |18j 
and [19] were not accounted for in |16| . As a re- 
sult, the values of form factors taken from [16] 
at the J/i/j and i/j(2S) are much smaller then the 
ones obtained in [17118119] and new investigations 
were necessary. 

For the pion form factor we keep the structure 
of the model from fl6l: 



F„(s) = 



N 



Y^c Pn BW n (s) 



.77=0 

oo 



J fit 



c Pn BW n (s) 

_n=N+l 



- 1 dual — QCD Nc = aa 

(5) 

where firsts N-l couplings are fitted (cn is cal- 
culated). We perform the fit taking into account 
the experimental data not accounted for in 16J. 
For kaons, unlike in [16j . we use infinite towers 
of resonances for p, uj and tj> (see [3D] for more 
details). The results of the fits are summarized 
in Fig|5j For pions we considered two versions 
of the model, where the Breit-Wigner function 
was taken at tree level [21] - Kiihn-SantaMaria 
(KS) model and where it was taken with pion loop 
corrections [32] - Gounaris-Sakurai (GS) model. 
For kaons we also considered two versions of the 
model, constrained and unconstrained (see |16] 
for the details). 

One has to remember that due to the finite 
detector resolution one never observes the true 
distribution of the events but its convolution 
with the detector resolution function. This is 
extremely important for studies of the narrow 
resonances, where typical energy resolution is 
much bigger than a width of a resonance. Fig[5] 
shows the differential cross section of the process 
e + e~ — > J/Y>7 — > 7r + 7r~7(7) with true distri- 
bution of the events (solid line) and the realis- 
tic differential cross section obtained by smearing 
with the Gaussian distribution with the 14.5 MeV 
standard deviation (taken from [2~4"]h 
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Figure 5. The experimental data |17|23j com- 
pared to the model fits results. The form factor 
at J/ip is from [TJ5] (theoretical extraction). 



We investigated also the role of the FSR radi- 
ation (at next to leading order) on the radiative 
return cross sections in the vicinity of the nar- 
row resonances. In Figl7]we show the results for 
two pion final state. Relatively big corrections to 
the differential Q 2 distribution, coming from FSR 
(IFSNLO=FSRNLO+ISRNLO) are observed as 
compared to the ISRNLO only, even if the inte- 
grated cross sections differ only by about 2%. 

4. Conclusions 

The present status of the PHOKHARA event 
generator was described. The four-pion channels 
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Figure 6. Differential cross section for y/s — 10.52 
GeV of the process e + e~ — > J/ip') — > 7r + 7r~7(7) 
without (solid line) and with the detector smear- 
ing effect. 



reanalysis and the latest upgrade of the genera- 
tor, implementation of narrow resonances, were 
presented. 
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